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We use bosonization methods to calculate the exact finite-temperature single-electron Green's 
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function of a spinful Luttinger liquid confined by open boundaries. The corresponding local spectral 
density is constructed and analyzed in detail. The interplay between boundary, finite-size and 
thermal effects are shown to dramatically influence the low-energy properties of the system. In 
particular, the well-known zero-temperature critical behavior in the bulk always crosses over to 
a boundary dominated regime in the vicinity of the Fermi level. Thermal fluctuations cause an 
■ enhanced depletion of spectral weight for small energies uj, with the spectral density scaling as lo 2 

for lu much less than the temperature. Consequences for photoemission experiments are discussed. 
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£j ■ I. INTRODUCTION 

c/3 ; 

In the last decade there has been an enormous interest in metallic phases of matter which are not Fermi liquids. 

S A paradigm for these is the Luttinger liquid, describing the low-energy, long-wave length limit of gapless electron 
systems in one dimension (lD)Era. The Luttinger liquid satisfies Luttinger's theorem, but the interaction wipes out 
the quasi-particle pole of the electron propagator, with the disjoint Fermi surface (consisting of two Fermi "points" 
±fci?) supporting only collective charge- and spin density excitations. These excitations are dynamically independent, 
effectively leading to a spatial separation of the charge and spin of an electron added to the. system. As a result, 
i ^ 1 1 spectral properties and dynamical correlations are quite different from those of a Fermi liquidLj. n n 

The generic low-energy properties of a gapless interacting ID electron system are well understoodtru in terms 
t-H , of universal power-laws of correlation functions apd spectral functions. The behavior of such so-called Luttinger 
liquids is coded in the Tomonaga-Luttinger model], much in the same way as the normal state of 3D interacting 
electrons is patterned on the free Fermi gas (Landau Fermi liquid theory). However, it is only in the last few years 
that laboratory technology has advanced to the point that the notion, of a Luttinger liquid can be confronted with 
experiments. Indeed, the ability to manufacture true ID quantum wires&El, as well as the development of high-precision 
_^ • spectroscopic techniques for probing quasi- ID materialsQ, have provided strong impetus for investigating Luttinger 
[~^. ' liquid physics in realistic contexts. Added motivation comes from the realization that the edge excitations of the 
0^ ■ fractional quantum Hall effect can be described in terms of a chiral Luttinger liquidQ. Also, some of the "non-Fermi 
liquid scenarios" for the normal state properties of the cuprate superconductors draw heavily from Luttinger liquid 
theory, suggesting possible extensions to higher dimensions!! 

Whereas,J-pttinger liquid theory has been successfully used to predict fractional quantum Hall edge state 
transportl!3liilll3, its applicability to the more traditional realm of quasi-lD materials remains controversial. This 
class of materials contains organic conductors such as tetrathiafulvalene-tetracyanoquinodimethanc (TTF-TCNQ) 
and the tetramethyltetraselenafulvalene (TMTSF^X Bechgaard salts with X = CIO4, PF6, etc. All available infor- 
£J . mation suggests thAt—the physics of these compounds is dominated by strong elactron correlations and pronounced 
one-dimensionalitycHla. Yet, data from NMRt3 and photoemission spectroscopy!] interpreted within a conventional 
Luttinger liquid framework seem to imply single-electron correlations governed by an exponent much larger than can 
be provided from any realistic lattice Hamiltonian (of which the Tomonaga-Luttinger model would be the effective 
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low-energy continuum theory). This poses a major problem for the modeling of these materials, and its resolution 
■ - - ■ remains a challenge for the theorist. 

One of the basic quantities to consider in this context is the single-particle spectral density, as this is the object 
that determines the outcome of a photoemission experiment. Surprisingly, thp^full calculation of the spectral den- 
sity of anjdeal, infinite volume Luttinger liquid was only recently performedcJ, expanding upon earlier results by 
SuzumuraEa. In the present paper we go a step further and investigate the local, finite-temperature spectral density of 
a Luttinger liquid confined by open boundaries (simulating a reflecting barrier or edge potential). As photoemission 
spectroscopy measurements are highly sensitive to boundary effects (with the photoelectrons traveling only a short 
distance, of the order of a few lattice spacings from the surface of the sample), it is crucial to incorporate an analysis 
of nontrivial boundary conditions. Also, spectroscopy on high-mobility quantum wires doped with artificial impurities 
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("antidots")£3 may soon be within experimental reach, adding yet another reason to study this problem: at sufficiently 
low temperatures a potential scatterer is expected to act essentially as a reflecting barrierEJ. 

In the present paper we calculate the exact finite-temperature single-electron Green's function and the associated 
local spectral density of a confined spinful Luttingcr liquid. The combined effects of boundaries and finite volume 
will be shown to strongly modify the well-known bulk Luttingcr liquid spectral density. In particular, the presence 
of boundaries causes a novel scaling behavior for energies close to the Fermi level, which produce a depletion of 
spectral weight significantly larger than that for a bulk system, as we have reported for a semi-infinite system at 
zero temperature beforelifl. In addition, thermal fluctuations deplete the levels further and give rise to a universal uj 2 
scaling of the spectral density for sufficiently small frequencies uj, which is analogous to the temperature effects in 
"bulk" systemaij. Particular attention is given to the nontrivial zero mode contribution to the correlation function 
from fluctuations in the spin and charge quantum numbers, which turn out to violate the spin-charge separation (i.e. 
the partition function cannot be split into separate spin and charge factors). Interestingly, the confined system is 
sensitive to the ratio between the effective velocities for the collective charge- and spin excitations (which in turn 
depend on the effective electron-electron coupling): As the velocity ratio locks into a rational value the spectrum 
separates into distinct peaks of equal spacing. The spacing between the peaks becomes dense for irrational values of 
the velocity ratios, leading to a quasi-continuous spectrum for large uj. This hints at a resonant interference between 
standing waves of charge and spin for special values of the electron-electron coupling. 

The paper is organized as follows: In the next section we introduce an extended version of the Tomonaga-Luttinger 
model and review its bosonization in the presence of open boundaries. In Sec. Ill the full finite-temperature Green's 
function in a bounded domain is derived, and the resulting local spectral density is extracted and analyzed in Sec. IV. 
Section V contains a brief discussion of possible consequences for photoemission experiments, as well as some con- 
cluding remarks. 



II. ID ELECTRONS IN THE PRESENCE OF OPEN BOUNDARIES: BOSONIZATION 

As our model we take an extension of the Tomonaga-Luttinger Hamiltoniarfl, which describes spinful Fcrmions in 
one-dimension with a (repulsive) local interaction. All electron-electron interactions must conserve spin and charge, 
so that we use the most general gapless Hamiltonian density that is invariant under the corresponding SU(2) and 
U(l) symmetries, 

W =V F ^L^^L^-^k^^R,' 

+ 92±JlJR a + 92\\JIJr + + W) + 9i\\{JlJl + JrJr) (1) 

where we have used the traditional "g-ology" scheme to index the couplingfl. The chiral Fermion currents are defined 
as ^l/r ~ tl/r a^L/R. <j ■> anc ^ ^L/R,a{ x ) ar e the left/right moving components of the electron field ^ a {x) expanded 
about the Fermi points ikp, 

9 v (x) = e-***1> Lta {x) + e^ R Jx). (2) 

This expansion is valid as long as the lattice spacing is much smaller than all length scales we want to consider. 
Therefore, the energy range around the Fermi-surface is limited to a region where a linear approximation of the 
spectrum is justified. This is often conveniently illustrated by a cut-off parameter, but more accurately one should 
take higher order operators in the Hamiltonian into account which result in corrections to the linear spectrum which 
are of higher orders in \k ~ kp\. These perturbations have non-universal coefficients which depend on the detailed 
interactions of the underlying lattice Hamiltonian. The appropriate value of the cut-off parameter (i.e. the range of 
validity) is then determined by the momentum scale at which those corrections become so large that the perturbation 
series no longer converges. Generically we can only roughly estimate the range of validity to be about one order of 
magnitude less than the band-width. 

The first term in Eq. ([!]) is that of free relativistic fermions, while gi and 174 describe forward electron-electron 
scattering. We have also explicitly included a backward scattering term gi±. The coupling constants depend on the 
microscopic parameters of the underlying lattice model and vf is the Fermi velocity. Normal ordering is carried out 
w.r.t. the filled Dirac sea, and we sum over repeated spin indices. 

Eq. (|l) defines a "standard model" for low-energy electrons in a ID metallic phase, and is easily derived from the 
Hubbard Hamiltonian 
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H 



H 



U>0. 



(3) 



In the weak coupling limit U <C i, we can treat f as a perturbation, and the tight-binding band e(k) — —2tcosak 
may be linearized around the Fermi points ±fep = ±n e 7r/2a, n e being the electron density and a the lattice spacing. 
The electron operators are replaced by the chiral fields ^ L j R a (x) in the continuum limit 



-ikpna 



As a result, the Hubbard model gets mapped onto the Hamiltonian density in Eq. (|I|) 
9i± = 92± = 2.94± = Ua and g 2 \\ = ff 4 y = 0. The Umklapp term e~ t4kFX tp[ a ^ R a ip f L _ -V'fl,- 



(4) 

with coupling constants 
+ /i.e. is also generated, 

but does not contribute away from half-filling (n e ^ 1, fcy/ ir/2a) due to rapid phase oscillations. For this case we 
are left with the theory in (jl|), with v F = 2at sin(fcpa). It is important to emphasize that the Hamiltonian in (|l|) 
faithfully represents the low-energy sector of the Hubbard model also for strong on-site repulsion iMH. However, when 
U is not small, the procedure above fails to identify the proper values of the model parameters, and instead these 
have to be inferred from the exact Bethe Ansatz solution of the Hubbard modelOE3. 

The Hamiltonian (Q) is conveniently bosonizeded by introducing charge and spin currents and the corresponding 
bosons 4> c and <j) s with conjugate momenta n c and n s , respectively 



J C r /S = —= [ Jl ± Jj ) = 



V2 



1 



47T 
1 



{dx<t>c/s - 
{dx<t>c/a 



n 



c/s) 



n 



c/s) 



(5a) 
(5b) 



The resulting theory describes separate spin and charge excitations moving with velocities (to lowest order in the 
coupling constants) 



#4|| 

v c = v F H 1 v s = v F 

7T 7T 

where v c > v s for repulsive interactions 174^ > 0. The Hamiltonian becomes 

U = {y[(^^) 2 + n '] + ^ [{d^uf-K]} + 5i± const, cos V8^0 S , 



9i± 

IT 



941 

n 



(G) 



(7) 



where g c = g 2 \\ + 921 and g s — g 2 \\ — g 2 ±- The charge interaction g c can be absorbed into the free Hamiltonian by a 
simple rescaling-of the charge boson, but the spin interactions g s and gi± obey Kosterlitz-Thouless renormalization 
group equationsEj with flow lines along hyperbolas gl — g\^_ = const, (to lowest order). For g s > —\g\±\ the spin 



sector develops a gap in the low energy, long wave-length limit, but for g, 
point g* a = — v/ffF 



< —\gi±\ the system flows to a stable fixed 
the interaction corresponds to one single marginally irrelevant 



9i± = 0- For 9s = -\9i 

operator, so that g* = g^ ± = 0. If the flow to a stable fixed point occurs, we can rescale the bosons by a canonical 
transformation to obtain a free theory [y = s, c) 



where to first order in the coupling constants 



K 2 a = 1 - 9 :/2irv s , K 2 C = 1 - . 9c /27TO c 



(8) 



(9) 



This yields the Hamiltonian 



(10) 



where v v = (v v + y~)K„, i.e. V v = v v to first order in the coupling constants so we will omit the 'bar' in the following. 

The chiral components of the electron field can now be expressed in terms of free boson fields and their duals by 
using Eq. (@) and the formulae^ 



1pL/R,A x ) 



_Jh_ 

\J2tux 



exp ^f2ni J J£^ R (x)dx 



(11) 
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where rj a obeys the anticommutation relation {?7j,?7x} = 0, with 77^ = 1. The presence of 77^ in (^Tj) guarantees that 
operators with different spin a obey anticommutation relations. 

Using the duality relation n„ = d Vv t<t>u — d x 4> Vl we thus obtain, using Eqs. (0) and (ph 



1pL/R,a(x) OC 7] a Y[ ex P ie ">°j^ (TKv4>v{x) - K v X 4>v{x) 



V — C.S 



(12) 



with e v a = 1 unless v = s and a =\ when it is equal to — 1 . Note that we have obtained the bosonization formula 
( |l2| ) with no assumption about boundary conditions. 

We now apply the formalism above to a system of length L with open boundary conditions and thus require the 
electron field \E' (T (a;) to vanish at x = and at x — L. From Eq. (j^) we see that this implies 

^i lff (0,t) = -^ R , a {Q,t) (13a) 
4>L,A L > t) = -e l2kFL ^ R AL, t) . (13b) 

Considering Eq. ( [l2|) this gives us fixed boundary conditions on the boson fields 4> v , which in turn determine the mode 
expansion. 

To calculate the mode expansion for the boson we find it most convenient to consider the classical Euler-Lagrange 
equation with fixed boundary conditions at x = 0, L 

cf> v (0,t) = C , cj> v {L,t) = C L , (14) 

and then performing a canonical quantization. We therefore consider the classical Lagrangian density corresponding 
to the Hamiltonian (10) 

£=Ey [(^>» 2 - [d^f] (15) 

V — S.C 

where 3 Vv t = ^-dt- The resulting Euler-Lagrange equations can be expressed in terms of d± = d x ± d Vv t<> 

{d 2 Vut -dl)<t> v = d + d^ v = Q, (16) 

and it follows that the two solutions can be written in terms of left- and right moving bosons, 4> v ,l(x + v vt) and 

4>vM.{. x — v v t). We use the combination <f> v (x, t) — 4>v,l{ x + v ft) + 4>v,r{ x — v v~k) an d its dual field 4> u (x, t) — 4>v,l{ x + 
v u t) — 4> u ,r{x — v v t), related by 

d x (j> v = d Vvt <f> v , d x <j> v = d Vt/t 4>u ■ (17) 

The classical solution with the boundary condition (Q) is obtained in a straightforward way, and after canonically 
quantizing we find the mode expansion for the quantum fields according to the boundary condition ( |l3| ) 

x sr^ 1 tvko. 
, — + > —^= sin — — 

n—l 



i/i / - X v-^ 1 . niTX ( - mvyut \ 

<Pv(x, t) — 4>v ,0 + Qu-r + > , — sin — — —le L o„+h.c. , (18a) 

L *■ — ' V717T L \ / 

n—l v 

00 

4> v {x,t) =4>„ fi + Q v -r + > ^=cos— — e _J ^^< + h.c. . (18b) 



The non-zero commutation relations among the mode operators are [<z^,a^'] = 1 and [<fiv,o,Qv] = i, while <j) U fi are 
c-numbers. This result (|l8| ) obtained from the classical solution automatically contains the correct zero modes, in 
particular the total charge and spin operators Q v , which were first postulated by Haldane for a periodic systemE. The 
resulting energy spectrum from Eq. ( |l0| ) is 



H =H I^ + E^n^ ' (19) 
v—s,c \ n—l / 

From Eqs. (jig) we can read off the mode expansions for left-moving bosons 
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<f> VjL (x } t) = - (j> v (x,t) +(f> v (x,i) 
= 2 (</v,o + <f>v,oj + Qv 



X + v v t 
2L 



£ 

n=l 



h.C. 



The right-moving boson field can be related to the left-moving one by 

<t> v ,R{x,t) = - 4> v {x,t) - 4> v (x,t) = -<j> U)L (-x,t) +4> v ,q. 



(20) 



(21) 



The boundary conditions ( |13| ) on the fermions provide us with the quantization condition for the eigenvalues of 
the operators Q v (using special commutation relations of <j) v and <p v at the boundaryc3 that follow from the mode 
expansion rtlq)) 



iir I 2k F L 

KcQc = \j^ U + 1 + — — 



K S Q S = x j —m , 



(22a) 
(22b) 



where m and n are either both odd or both even integers. Therefore, the quantum numbers for the total spin m and 
the total charge n are not independent, which is to be expected because we can only insert and remove real electrons, 
i.e. we cannot change the total charge and the total magnetization independently. In this sense these degrees of 
freedom do not obey spin-charge separation and the partition function does not factorize. It is interesting to note 
the formal similarity to a spinless periodic system, which is also described by two channels (left and right moving) 
and where a similar condition holds for the total current and charge quantum numbersa. As we will see later, the 
zero modes indeed give a contribution to the Green's function which does not factorize, while all dynamical degrees 
of freedom in Eq. ( |l9|) remain spin-charge separated. 

The constants Cj o = \p^K~ 1 and (j) s .o — are also determined by the boundary condition (modulo the intrinsic 
periodicity of the boson \[2~kK~ 1 ) . Hence, we obtain from Eqs. (13) and (pl|) 



1pR,a(x,t) = -1p Lja .(-X,t) . 



(23) 



which allows us to write the full theory in terms of left-movers only. This concludes our analysis of bosonization in 
the presence of open boundaries. For an alternative approach - exploiting the path integral formulation - see Ref. f27L 



III. GREEN'S FUNCTIONS 



Using the formalism above, the exact single electron Green's function for a confined Luttinger liquid with open 
boundaries can now be calculated. With the decomposition in Eq. (|^) and using Eq. (p3|) , we have 

<*£ (z, t )9 a (y, 0)) = e ik '^G{x, y, t) + e -^ (x '^G(-x, -y, t) 

-e ikF{x+y) G{x,-y,t) ~ e- ikFix+y) G(-x,y,t), (24) 

where the chiral Green's function G(x, y, t) = (ip^ a (x, t)ip L a (y, 0)^ is derived in Appendix A. The result is a product 
of the spin and charge contributions F s c and a factor H from the zero modes 



G{x,y,t) oc H(x,y,t) ] ] (F v (v v t + x - y)) ' " ' (F u (v v t - x + y)) ^ ' 



(25) 



|F„(2aO||F„(2y)| 



F v {v v t + x + y)F v {v v t - x - y) 

The contribution H from the zero modes is given by (see Appendix A) 

H( i \ _ + T c k F L\T c )'d 3 (u s \T s ) + i} 3 (u c + T c k F L\T c )t9 2 (u s \T s ) l2Uc h£± 

V'W- Mrck F L\r c )M0\r s )+Mrck F L\r c )M0\rs) 



(26) 
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where the Theta functions, i?fe(u|r), are defined in sections 8.18-8.19 in@ and 



7T (v v K v H + x- y\ ,v v fi 



2 V L ) K*L 

The expression for,. H does not contain any poles, but may still influence the spectral properties significantly in 
mesoscopic systemsotj. It is also interestingly to note that the contribution from the zero modes cannot be written 
as a product of independent spin and charge contributions, which is expected since the quantum numbers n and m 
in Eq. ( p2f ) are not independent. 

The factors F v frpap the spin and charge bosonic modes can also be written in terms of the Theta functions (see 
sections 8.18-8.19 i: " 



F v {z) = i— sin^ TT 



an 2L 

fe=l 



sin ff 

sinh k ™L 



— n2Ll 2L i . (28) 



Here x and y denote the distance from the left boundary (a; = 0) and the argument z carries an implicit cut-off z — ia. 
The parameters v v and K v are defined in (^J) and @ respectively. 

We can immediately verify that the full Green's function is anti-periodic under translation by the inverse temperature 
t — > t + i/3 as it should be. It is also interesting to note that the factor from the zero modes is periodic under the 
change of the Fermi-level hp by ir/L, which results in periodic oscillations of the spectral properties of the system 
as the chemical potential is changed. This is a_tnanifestation of Coulomb-blockade oscillations, i.e. resonances can be 
observed if the system is of "mesoscopic" sizeEJ. It is important to emphasize that the zero mode contribution H in 
Eq. ( prf ) is obtained by using a grand canonical ensemble when taking the averages, thus allowing for fluctuations in 
the magnetization m and in the particle number n. By fixing these quantum numbers (i.e. using an idealized "closed" 
system) or by letting the system size tend to infinity, the zero mode contribution collapses to a constant phase. This 
provides a vivid example of how different statistical ensembles may lead to different results on mesoscopic scales and 
downwards, where quantum coherence effects become important!!! 

We see in Eq. (|2q ) that we recover the universal power-laws, which give the expected branch cuts in the Green's 
functional. In addition, we get a contribution from the boundary, which is entirely contained in the last factor of 
Eq. (^5|) and gives an additional analytic structure. This factor does not contribute in non-interacting systems 
(K c = K s = 1 ) , so that in this case the presence of the boundary is seen only by the addition of the two last "Friedel" 
terms in Eq. ( |24| ) compared to the bulk case. In contrast, with interaction (K s , K c < 1), the boundary influences also 
the chiral Green's function. This is expected, since only with interactions can effects from electron scattering off the 
boundary propagate to other parts of the system, thus influencing also the chiral pieces of the full Green's function. 

The zero-temperature limit T — > is readily obtained by letting j3 — + oo in Eqs. (|2^) and ( p8|) 



G{x,y,t) oc e 4(2n °- 1)Uc ]Q — sin- 



2L . n(v v t + x — y) \ ( 2L . ir{v v t — x + y) 



TT 2L I V 7T 2L 



sin ■ 



sin H2L S i n m 



■K(x-\-y-\-v v t) • Tr(x-\-y—v u t) 
Sill 2L sm ?L 



(29) 



where no = {kpL/i: mod 1) effectively measures the difference between the Fermi- vector and the highest occupied 
level (which are not necessarily the same in a system with discrete energy levels). The phase e l ( 2Tl o-i)«c comes 
from the zero modes and does influence the time correlations. For special values of no this phase may have a different 
dependence on u v . When this phase is neglected, Eq. ( p9|) agrees with the resultsjn Ref. obtained via bosonization 
according to the "Haldane prescription"!!!. It is also in agreement with resultsta, obtained by conformally mapping 
the semi-infinite complex plane onto a finite strip. 

Using the Poisson summation formula, we can obtain the limits of the Theta functions i9i(z|i7) w 
2 7 -i/2 e -z 2 A7 e -T/4 7 sinh * and ti 2 {z\i"f) « #3(^7) ~ 7 -V2 e -2 2 A7 as 7 0. Hence, by letting L -> 00 in Eqs. @ 
and (|28|), we obtain the finite temperature chiral Green's function of a semi- infinite system with open boundary 

n , ,x tt (v v P . ,n(v v t + x - y) \ 5 (v v f3 . ir{v v t ~ x + y) " 
G(x, y, t) oc J] ( — -nh — \ I— smh -g— 



G 



Finally, the T — ► limit for a semi-infinite system is obtained by letting (3 — > oo in Eq. (30) (or L — > oo in Eq. (p9|) ) : 



G(x,y,t) oc J] M + ar-y)" • (M-x + y)" - ( [(iE + y /_ ^ j ( 31 ) 

We note that in the limit xy ^> \(x — y) 2 — u£t 2 |.-jthe last factor in equation goes to unity and we recover the 
known zero-temperature bulk correlation function^ (as we also do in the non- interacting case K c = K s = 1). In the 
limit of equal time, t = 0, this gives the asymptotic correlator 



1 

(x - y) 2 



G(x,y,t = 0)<x 7Z A (32) 



with the equal-time bulk exponent A = (K 2 + K 2 + K c 2 + K s 2 )/8. In contrast, when one of the points is close to 
the boundary compared to the relative distance x 3> y one finds 



1 

(x - y) 2 



G(x,y,t = 0)<x - (33) 



with Aj_ = (3K C 2 + 3K S 2 + K 2 + K 2 )/16. This is to be compared to the asymptotic long-time behavior of the 
autocorrelation function in the limit t 3> x,y, which behaves as 

G(x,x,t) oc (34) 

with A|| = (K~ 2 + K~ 2 )/4. As expected from scale invariance, we thus recover the scaling law A^ = (A + A||)/2. 
It may be worth pointing out that the dynamic (t ^ 0) asymptotic large distance correlator in Eq. (|33| ) is not 
governed biz_the single exponent Aj_, as one may naively have expected from the analogy with classical critical 
phenomenaEj. Instead it remains a product of separate charge- and spin correlators with exponents (3JsT~ 2 + K 2 )/16 
and {3K~ 2 + K 2 )/16, respectively. This is not in conflict with scale invariance since the theor y is built out of two 
distinct sectors, each with its own effective velocity. More importantly, the behavior in Eq. (gij) reveals that the 
asymptotic low-energy behavior of a Luttinger liquid with an open boundary belongs to a different universality class 
than that of the bulk theory. We shall elaborate on this in the next section when we discuss the local density of states. 



IV. LOCAL SPECTRAL DENSITY 



To understand the physical implications of the boundary correlations we study the local spectral density N(u,r), 
given in terms of the single electron Green's function in ( gjj ) 

1 f°° 

N{u,r) = —J e^ f (n 0), * CT (r, t)}> dt, (35) 

where uj is measured relative to the Fermi energy and r is the distance from the boundary. 

At T — and without the boundary, the integral in ( |3E| ) can be done exactlyfij and one finds that the spectral 
density scales at the Fermi level as N(uj) oc uj a baik^ w here the exponent in the bulk is given by 

«buik - (Kc + K- 2 + K 2 + K; 2 )/A - 1. (36) 

However, the boundary clearly influences this scaling behavior, and by inspection of Eq. ( |3l| ) for a semi-infinite 
system at T — 0, simple power counting reveals that there must be a crossover to a boundary dominated regime for 
ruj < v ci v s with a novel exponent 

abound = {K- 2 + K- 2 )j2 - I . (37) 

Interestingly, the boundary exponent abound therefore always dominates for sufficiently small u>. It is also interesting to 
note that the last two terms in equation (E4J) make a contribution which oscillates at twice the Fermi wave- vector and 



7 



drops off with the distance from the boundary proportional to e l2k Fr r -(K a 2 +k s 2 )/2^ ^his contribution is reminiscent 
of a Friedel oscillation, although it can probably not be observed directly, since experimental measurements of the 
density of states (in particular photoemission) effectively averages over several lattice sites. We therefore ignore those 
"Friedel" terms in the following calculations and make the replacement 

<{*t ( r , 0), * ff (r, t)}> G(r, r, -t) + G(-r, -r, -t) + G(r, r, t) + G(-r, -r, t) , (38) 

in Eq. (|3^). Using our exact results for G(x,y,t) in the previous section allows us to fully explore the physically 
relevant piece of the local spectral density. 



A. The zero temperature and infinite length limit 



For T — > and a semi-infinite system, L — > oo, it is readily derived from (plh, (pa) and 



that 



N{uj 7 T = 0,r) 
where 



rft cos7(t) 



coswi — 1 



2r 



-bo/2 



Vat 

2r 



-b B /2 



(39) 



§(a s + a c ) 
7 (t) = { § (a, + a c + 6 c ) 

f (a s + a c + 6 S + b c 



< i < 2r 

— < t < OO 



and 



Ki + K- 



K = 



K- 



K 2 



(40) 



(41) 



4 4 

To cure the original divergences in (p9|) we have subtracted an infinite constant by renormalizing the static spectral 
density to zero for any given r, i.e. N(lo = 0, T = 0, r) = 0. 

Let us study the two limiting cases r — > a (boundary regime, with a the short-distance cutoff) and r — > oo (fettZA; 
regime). Both limits are easily obtained from Eq. (p3) by using the integral (3.823 inc3) 



dx(cos x — l)x 



-k 



1 



T(fc) 2cosf k 



1 < k < 3. 



with r(fc) the Gamma function. This gives 



N{oj) 



7iT(l+a b ulk) 



7rr(l+Q b 







T — > OO 

r — » a 



(42) 



(43) 



where oibuik = a s + a c — I and pbound — a s + i c + & s + b c — 1 . In the bulk limit r — > oo this result is in full agreement 
with previous calculationsEJ'E3, but in the boundary limit r — > a we observe completely different exponents. By 
substituting uxt by x in the integral in Eq. ( |39| ) we see that apart from a prefactor measuring the distance to the 
boundary, the spectral density depends only on the scaling variable ru). This implies that the condition for boundary 
behavior is ^ -C 1, i.e. regardless of the value of r there will always be a region in lu around the Fermi-energy where 
the spectral density is determined by the boundary exponent abound- Actually, the scaling behavior of the spectral 
density splits into three distinct regions in lu where different exponents govern the leading scaling. 

As an example, let us choose parameters adapted to a description of the large-[7 Hubbard chain away from half- 
filling (cf. Sec. II and Ref. |^). In this case-tte SU(2)-invariance forces K s = 1, and it is known from Bethe ansatz 
calculations that — > 1/2 when U — > oocllfi3. Therefore, the spin channel is not affected by the boundary (since 
b s = 0), and the local spectral density N(u>,r) in Eq. (|39|) splits into two asymptotic sectors only: The boundary 



regime for u> <§C ^ and the bulk regime for u> ^> — . From these numbers, the well known result N(ui) oc 



follows 



immediately for the bulk regime, as can be seen from Eq. (|43j). In the presence of the boundary, however, we cross over 
to the boundary exponent abound = 1/2 for to < v c /r. The results in Fig. (|l|) clearly show the crossover from boundary 
behavior for ru>/v c < 1 with exponent abound = 1/2 to bulk behavior for rco/v c > 1 with exponent abuik = 1/8 (the 
corresponding power-laws are superimposed). In the figure, the distance from the boundary, ?*, is held constant, thus 
setting the scale. The observed oscillations in Fig. (Q) are an intriguing secondary effect, which vanish asymptotically 
as sin(2wr/u c )((jjr) bc / 2_1 . It is important to emphasize that they are not due to the "Friedel" terms in equation (|24|) 
which have been neglected. Instead, they originate from the integrable singularity of the integrand in Eq. (|3^) at 
t = — , which is only present in the boundary case. 
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B. The finite temperature and infinite length limit 



We now consider the effect of finite temperatures on a semi-infinite system with a boundary. In general, by turning 
on temperature one induces a different behavior for small ui: Both in bulk and boundary regimes (defined as above) 
the spectral density crosses over to uj 2 scaling when uj < 2w//3. This is due to thermal fluctuations which produce 
exponential damping of the density correlations for unequal times. The formal expression for N (uj, /3, r) for this case, 
with 3 — 1/ksT and fee = 1 is given by 



N(<jj,/3,r) = — xv~ a "v~ a ' / dt cos-y(t) 



cos uit 



/smh^t 



sinh ^ (2r + v s t) sinh ^ (2r - v s t) 



sinh 



2 2ti 



-b s /2 



sinh ^ (2r + v c t) sinh ^ ( 2r - Vc t) 



sinh 



2 2ti 



2r 



-6 c /2 



-bc/2 



2r 



-f>„/2' 



(44) 



where 7(f) is given in Eq. (40). We have subtracted the same infinite constant as in Eq. (|39j), thereby, as before, 
renormalizing the zero-temperature static spectral density to zero for any given r, i.e. N(u) = 0, T — 0, r) = 0. 

Substituting uit by x in the integration of Eq. (44) as above reveals that N(ui,8,r) can be expressed as a function 
of two scaling variables rui and uid (up to the same r— dependent prefactor as for zero temperature). By inspection 
we recover the T — result when ^ 3> 1, but for small ui a new behavior sets in. Independently of whether we are 
in the boundary or bulk region, N(u), (3, r) ~ N(0, f3, r) will be proportional to ui 2 for small u>(3 



N(u, (3, r) — N(0, f3, r) oc / dt(coscjt - 1) sinh -t 

Jo VP 

3 f°° . k . 2^0 wjSmmU 

— —2— / ax(smhx) sm — x — > 

TT./n ^ ^ 27T 



2 1 //3 
-a; - — 

2 \ 7T 



da: (sinh x) k x 2 (45) 



where the last integral converges if < k < 3 (in our case k = a n + 1 where a„ is the boundary or bulk exponent). 
In conclusion, the spectral properties are unaffected for energies well above the temperature ui ^> =£■, as expected. 

However, the spectral density will exhibit a uj 2 behavior for uj < ^S- before the cross-over to the T — behavior occurs, 
as can be seen in Fig. (|^) where we again have considered the large-U Hubbard model away from half-filling. This is 
in complete agreement with the recent work by Nakamura and Suzumuralij where the analogous effect was reported 
for an infinite "bulk" system. Effectively, Fig. (^) contains all information about both the bulk and the boundary 
case, since we are free to adjust the distance from the boundary r to any value and this only changes the scale on 
which we measure the energies and temperatures. Therefore, we observe a cross-over from the quadratic behavior 
directly to bulk behavior if |p > 1, while an intermediate boundary region can be observed for ^ < 1. As we can 
see in Fig. (^), the spectral density can look very flat around the Fermi-level in either case, and the sharp cusp which 
has been predicted for T=0 may not at all be visible in experiments. 

It is interesting to note that the boundary exponent also shows up in the temperature dependence of the static 
spectral density N(u> — 0,3,r). This is expected since the static density samples all times, with the asymptotic 
large-time behavior governed by the boundary exponent. As can be obtained from (^J), 



N(uj = O,0,r) 



r « h o» nd (^)6 c (^)b sC .( abound + 1} 



— «: 1 

^> 1 



(46) 



where 



C(k) = cos — k 



dx ( sinh 



-k -k 

x — x 



(47) 



which is convergent for 1 < k < 3. With parameters again chosen to describe the large-£7 Hubbard chain (away from 
half- filling), the boundary dominated regime opens up for < |g < 1, as depicted in Fig. (§. 
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C. The zero temperature and finite length limit 



We now turn to a confined system with open boundaries at both ends. At this point we want to emp hasize that 
there is a distinction between effects that arise from a non-trivial boundary condition (as discussed in Sec. [V A) and 
effects from a finite system size (which may or may not have trivial boundary conditions). 

A confined system with open boundary conditions is technically more difficult to analyze, since the function F v (z) 
in Eq. (28) is periodic in z with twice the system size 2L. Because of the non- integer exponents in Eq. (^9|), we need 
to carefully keep track of the overall phase as we integrate around the various branch-cuts. However, we can simplify 
things by explicitly using the short distance cutoff z — > z — ia, which allows us to make a Taylor expansion of the 
factors in Eq. (E9), 



?' 2L£ \ 7 oo 
7TZ\ 7 / e 2' 

I sin — I 



/ 7TZ 

V S[n 2LJ 



(48) 



with 



cud) = (-iy 



r(7 + i) 



T(n- 7 ) 



r(n + l)r(7 - n + 1) r(-7)r(n + l) 



(49) 



An immediate consequence of the periodicity is that to gets discretized, which is consistent with the appearance of 
discrete energy levels for a finite system. The integral over the exponentials in the expansion (^8|) will give delta- 
functions at those special values of ui and we can try to extract an effective behavior in the prefactors, i.e. the 
coefficients c n . For a single channel case we can verify by inspection that the asymptotic behavior of the prefactor in 
Eq. ( f49| ) gives the expected power-law for large n, i.e. the semi-infinite length result can be recovered. 

However, for two channels and arbitrary values of lu we need to make a more careful analysis. By using the 
multiplication formula 



^2 a kX k ■ ^2 bkxk 

k=0 k=0 



E 

fc=0 



c k x 



Cn — ^ OfAi-fc 



(50) 



fc=0 



we derive from (\2 



§), ©, and m 



N(u,r) = — l—) 
air \ L I 



a s +a c 



2 sm — 
L 



b s +b c 



EE 

n=0 m=0 



f n (a c ,b c ,r)f m (a s ,b s ,r) 



nv c (a c + b c ) 7TV S (a s +b s ) 



2L 



2L 

TTV c {a c 



2L 



1TV C TTV< 

n m — - 

L L 

irv s (a s + b s ) 



2L 



(51) 



TTVr 



where 



f n (ki,k 2 ,r) 



p=0 



P 

-E 

9=0 



E c q{-k 2 /2)c p - q (-k 2 /2) cos 



„ 2nr 
(p-2q) — 



(52) 



To understand the role of the delta- functions in ( |51| ) it is convenient to represent the argument of the delta function 
in a two-dimensional parameter space, coordinatized by the pair of summation indices (n, m) in j5l|), as shown in 
Fig. (^J). The line connecting (lu, 0) and (0,o;) indicates where the argument of the delta- function vanishes, and hence 
selects the terms to be included in the double sum in (|5l|). The points are the allowed values of uj and the ratio v s /v c 
determines the relative distances between the points in the x- and y-directions. This ratio plays a crucial role, because 
when v s /v c is a rational number we have a resonance situation and the spectra will consist of peaks with constant 
spacing. On the other hand, for v s /v c irrational, N(u>,r) is still discrete for small lo but approaches a continuous 
function for large u>, since the number of points close to the line increases with increasing to. Moreover, if the spin- wave 
velocity v s is significantly smaller than v c the spectrum may appear continuous, but the discrete charge peaks may 
still be resolvable. (As v s — > the charge excitations are described by yet another exponents.) 

When the (experimental) energy resolution Aw is larger than the spacing between the peaks, it is appropriate to 
convert the infinite sums in Eq. (Bfi) into integrations over continuous variables. The resulting double integrals can 



10 



actually be done exactly in the extreme boundary case as well as in the extreme bulk case, and the semi-infinite length 
results are recovered in both scenarios. Even for the intermediate case it appears that the coefficients ( |49|) under the 
double sum closely reproduce the power-laws of the semi- infinite case if some smearing is taking into account. This 
is in strong contrast to the exponents that describe the jaomentum distribution, which are known to be strongly 
influenced by finite size effects as well as boundary effectslla. 

In conclusion, the main difference between a semi-infinite and a finite system is therefore the appearance of a 
discretized spectrum, with possible beatings of charge and spin excitations. This effect can only be observed in very 
small (mesoscopic) systems or with a very high experimental resolution, since the smeared spectral weight appears to 
follow the same frequency dependence. The effect of the boundary remains dominant for small frequencies in either 
case. 



D. The finite temperature and finite length limit 



The most general case is to consider both finite temperature and a confined system (finite length) . The periodicity 
of the Green's function (|2^) is unchanged, but we expect that the coeffic ients in front of the delta-functions will acquire 
temperature dependent corrections similar to the ones discussed in Sec. IV B. We can make a similar Taylor-expansion 
as in the previous section by using Eq. ( [so]) to expand the temperature dependent factor in Eq. (|2^) to the power of 7 



n 

fe=i 



where 



2L 



sinh k 



(53) 



fc=l 71=0 



d„( 7 ,z) = (-l)»^ c „_ p (-2 7 )^ Cg ( 7 )c p _ g ( 7 )e*-^ 2 

p=0 q-0 



and the coefficients c n are defined in Eq. (|49|). This can be written as 



(54) 



n 

k=l 



2L 



sinh 



oc m 
m=l n=0 



-m0- 



(55) 



where g n ,m{l) is a highly nontrivial but well-behaved function composed of the coefficients e n ( 7 ). The zero mode 
term H also makes a non-trivial contribution which will be discussed elsewhereEj. 

The integration over the exponentials in Eq. ( |55| ) gives us again delta functions. However, there is no shift or 
smearing of the peaks in the spectra due to temperature, only the height of the existing peaks are modified, i.e. the 
points where the delta- functions contribute are at the same values of lo as indicated in Fig. (Ji|). In the extreme 
boundary case we can derive an explicit (but complicated) expression for the spectral density and we observe that the 
temperature has a negligible effect for large w as expected. We conjecture quadratic behavior of the coefficients for 
small lo, which is supported by preliminary numerical evidence, i.e. we observe a similar effect to the one discussed in 
section IV B. Thus, we are left with the analogous conclusion from the previous section that finite size effects always 
result in a discrete level spacing of delta-functions, but do not alter the (smeared) dependence on freq uency . We 



therefore recover the same cross-over from co -behavior to boundary or bulk behavior as discussed in Sec. IV B. 



V. DISCUSSION 



In conclusion, we have derived an exact closed-form expression for the single-electron Green's function of a spinful 
Luttingcr liquid at finite temperature and confined to a finite interval by open boundaries. By analyzing the cor- 
responding spectral density we have obtained detailed information about the interplay between boundary, finite-size 
and temperature effects in an interacting electron system. Most importantly, we find that the scaling of the zero- 
temperature spectral density with frequency lo close to the Fermi level is always governed by a coupling-dependent 
boundary exponent significantly larger than the bulk exponent. In other words, the asymptotic low-energy behavior 
of a Luttinger liquid with an open boundary belongs to a different universality class than that of the bulk theory. 
Thermal fluctuations at finite temperature T > destroy this behavior and open up a new regime for lo less than the 
temperature where the spectral density exhibits quadratic scaling in lo. Not surprisingly, the same effect is present in 
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a bulk system for this frequency rang , implying that the boundary plays no decisive role in the process. In the case 
of a finite interval confined by two open boundaries, our results reveal a discretized spectrum with delta-functions 
at the allowed energy levels. It is interesting to observe the dependence of spacing of the energy levels on the ratio 
v c /v s between the effective velocities v c and v s of the charge- and spin excitations, respectively: For v c /v s a rational 
number, the spectrum consists of well-separated peaks which, for sufficiently large ui, coalesce to a quasi-continuum 
if v c /v s is shifted to an irrational number. Although this effect can not be observed experimentally, it nonetheless 
suggests a resonance phenomenon with the collective charge and spin excitations showing interference effects at special 
values of the electron-electron coupling. However, despite the appearance of a discrete spectrum we find that finite 
size does not influence the lu dependence of the integrated (i.e. smeared) spectral density significantly, so that the 
same boundary and finite temperature effects as for a semi-infinite can be observed. 

Let us close by briefly discussing the possible relevance of our results to experiments, in particular the photoemission 
studies on the Bechgaard saltsQ referred to in the introduction. These materials are composed of molecular chains 
which become conducting above some characteristic temperature T c , and are expected to show Luttinger liquid 
behavior provided the temperature is high enough to mask the weak inter-chain coupling. However, as mentioned in 
the introduction, high-precision photoemission experiments indicate a scaling of spectral weight with frequency that is 
inconsistent with standard theory of a bulk Luttinger liquid: The effective exponent a for scaling of the-photoemission 
intensity is roughly 1.25 (Ref. 0) as also seen in independent NMR experiments on the same materials!!-], whereas the 
largest realistic, value obtainable from a bulk Luttinger liquid description is a = 0.125, corresponding to_the large-U 
Hubbard chairjj. Attempts to include long-range Coulomb-repulsion, which can be shown to increase c£j, fails due 
to the instability against an insulating phase at a = 9/16E3, so other explanations must be invoked. As the typical 
escape depth of photoelectrons in the UV range is only 5 - 10 A, the experiments are extremely surface sensitive, 
suggesting that ID boundary effects may play a role for the observed scaling behavior. 

Consider first the case where one probes electrons that escape from a crystal face perpendicular to the ID molecular 
chains. The photoemission intensity I{lu,[3) is then proportional to the local spectral density N(ui, f3,r), integrated 
over the escape depth of the photoelectrons, and weighted by the Fermi-Dirac distribution fpniftoj) 

I(w,P) oc J dr f FD (pLu)N(uj,{3,r) . (56) 

In a boundary dominated region, I(ui) is seen to be dramatically reduced compared to the a bulk regime, considering 
our results in Eq. (|l3|). With a typical escape depth of a few lattice spacings, the condition for boundary behavior 
TLojv c < 1 may apply over an energy range of several hundr.ed meV (since v c > aEp, with a the lattice spacing and 
with Ep w 0.5 — 1 eV, depending on the particular materialEj). 

In the recent-photoemission experiments on (TMTSF^PFe (Ref. |7|), the chains are always in the plane of the 
cleaved surfaceO, and it is less clear to what extent ID boundary effects contribute. However, in the likely case that 
the cleaving of the surface introduces defects in neighboring chains, effectively breaking these into smaller segments, 
we may model the breaks by open boundaries and apply our results. Unfortunately, the actual defect concentration 
remains unknown, and it is therefore difficult to make a quantitative prediction from our results. This is an important 
issue that in principle should be possible to resolve via STM techniques. To explore the size of the boundary effects 
experimentally, it would also be of great interest to do photoemission experiments on cleaved surfaces that are 
perpendicular to the chains that could then be compared to the results from cleaved surfaces parallel to the chains. 

As discussed ir£3, the finite-energy resolution of the photon lines effectively introduces an averaging over the "true" 
spectrum 

I(u) ohB ee — L- / e-< w - ) 2 / 2A2 /(x) dx (57) 

which completely wipes out the power-law singularities in either the bulk- or boundary case. This "smearing" results 
in a similar effect as the thermal fluctuations, which also wipe out the sharp cusp from the power-laws as shown 
in Fig. (H), but have to be taken into account separately. With an experimental resolution of A = 20meV and 
at a temperature T — 50K (experimental values according toQ), and assuming boundary dominated behavior for 
N(lu, /3,r), the observed intensity in the vicinity of the Fermi level indeed appears to be depleted with an exponent 
of one or larger as shown in Fig. (||). (This should be compared to the large-U Hubbard exponent abuik = 1/8 
of the bulk spectral function without temperature or averaging effects). In experiments the condition for boundary 
behavior lu < v c /r will be satisfied over an energy range lu ~ Ep/r around the Fermi-level, where r is the distance 
from the boundary in units of the lattice spacing. This means that if the broken chains close to the cleaved surface 
have an average impurity density of a few percent, boundary effects could be observed over a region of up to lOOmeU 
around the Fermi energy. Experiments indeed suggest a scaling I(u>) b s oc oj a ° bs with a b s > 1, extending, howe\f«|, 
over a larger energy range. Thus, some additional mechanism (inter-chain coupling or electron-phonon couplings) 
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most likely have to be invoked to fully explain the data. Yet, a complete modeling of photoelectron spectroscopy on 
quasi-lD organic metals must certainly incorporate the boundary and temperature effects predicted in the present 
paper. 
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APPENDIX A: CALCULATION OF THE CHIRAL GREEN'S FUNCTION 

Here we extend the calculations of Ref. ^ to the case of spinful Fermions with open boundaries. To calculate the 
chiral fermionic Green's function, we find it useful to treat the contributions from the zero modes and the dynamic 
bosonic modes separately. Hence we write <j) u ,L in Eq. ( po|) as a sum of the zero modes and the harmonic oscillator 
terms 



0u,0 



X + v u t 



S u 



where the bosonic operators are contained in the sum 

1 



E 



■ n7r(x + u„t) 

e- 1 



(Al) 



(A2) 



We now insert this mode expansion into the bosonization formula ( |12| ) for ?pL, and use the definition of G(x,y,t) 
<>,*)V>i, CT (y,0)) to find 



-— (v v t + x-y) 



G(x,y,t) ocH(x,y,t) exp 

'-(K- 2 - Kl) (2B VtL {x, t- -y, 0) + 2B VjL (-x, t; y, 0) 



exp 



~{K v -KZ l ) 2 B v , L {-x,t;-y,0) 



'•■••V \ -{ K » + K-^fBvxix^t-y,®) 

4 

B ViL (x, t] -x, t) - B v , L {-x, t; x, t) - B VtL (y, 0; -y, 0) - B u , L (-y, 0; y, 0) 



x exp 



(A3) 



Here we used the identity e e 
bosonic Green's function 



-: e 



AB+- 



for the bosonic operators S v l and we have defined the 
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B^ L (x,t;x ,t ) = {S v ^ L (x,t)S v ^ L {x ,t ) - - [S^ L (x,t)S^ L (x,t) + S u ,l{x ,t )S v> l{x ,t )] 
The contribution from the zero modes is 

H(x,y,t) 



(A4) 



(A5) 



where u,. 



2L 



(v w K~ 2 t + x — yj. This factor cannot be written as a product of spin and charge expectation values 
separately, because the quantum numbers n and m in Eq. (|2^) are connected by the condition that both are even 
or both are odd. However, since we know the quantization condition ( |2^ ) and the energy spectrum (^|) for the zero 
modes, we can directly sum over all eigenvalues m, n. The factor H can then be expressed in terms of the elliptic 
theta functions (see sections 8.18-8.19 ino) 



„, ,s _ $2(U C + T c k F L\T c )l) 3 (u s \T s ) + ^ 3 (m c + T c k F L\T c )-d 2 (Us\T s ) i2Uc h£± 

^ y >>- d 2 (T c k F L\T c )§ 3 (0\T s )+MT c k F L\TM0\r s ) 6 ' " ' 



(A6) 
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where t v = i— - - 



To calculate the bosonic Green's function -B„.l, we insert the expression (A2) for S v ^{x,t) into Eq. 
gives 



1 

Awn 



e -2mn-i ^ >- 



1 (1 + m v n ) + ^ 



1 ml 



(|A4|), which 
(A7) 



where are the Bose-Einstcin distributions 



in' = ( a^a 1 " 



e' 3 ^ - 1 



ri -Wjyt + a; — t +a; ) 0*0%/ "K 

We define a = e~ 5l and 6 = e p_ ^, which allows us to write 



(A8) 



1 



1 -(&-*)' 



$>-")* = 5>- fc y 



fc=0 



fc=0 



The bosonic Green's function can then be written as 



oo 

^ 4nn 



n=l 
oo oo 



(a»-l)£(&-T+(«- B -l)r B £( 6 " fc )' 



fc=0 



fc=0 



fc=0 n=l 



oo oo 



We now can use the formula 



k=l n=l 



oo t. 



47T7J 



fe=i 



and we need to use the high momentum cut-off a in the first two terms of Eq. (A10) 



OO OO 

^_( a »_l) — (a"-l)c n , c = e~ Q £, a^O. 

' * 71 ' * 71 



This yields 



71=1 



B„, L =-— In 



4-7T 



47T 



1-acp (l-a6- fc )(l-a- 1 6- fe ) 

A' 



i - c AJ- (i - b- k )(i - b- k ) 

(ac) 1 ^ f ( aC )-l/2 _ ( ac )V2 ^ ( x + ( 2z j 



2L 



2?: 



n 

fc=l 



b k/2_ b -h/: 



Inserting a and 6 defined above gives us 



u I + i + '\ i ( vJ + x-vJ' -x' \ 1 
B v , L {x,t;x ,t ) = - I )-^ ln 



where 



2Z/ 7TZ TT 

F„(z) = i — sin - |[ 
CK7T 2L - L - L 



fc=l 



sing 



sinh fc 



TTV u f3 



F v {v v t + x — v v t' — x') 



The first term in ( A14 ) cancels with the phase in the zero mode part of (AS), resulting in Eq. (J29) 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 
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FIG. 1. The spectral density as a function of ruj in arbitrary units (from Ref. |l8[ ). The corresponding power-laws for 
«t>uik = 1/8 and abound = 1/2 are also shown. The distance from the boundary r is held constant and just fixes the scale. 
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FIG. 2. The spectral density as a function of ru> in arbitrary units for different temperatures r//3. For u < 1//3 we have 
parabolic behavior which crosses over to the T = behavior for u> 3> 1//3. 
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FIG. 3. The spectral density at uj — as a function of i n arbitrary units. The corresponding power-laws for abuik = 1/8 
and abound = 1/2 are also shown. The distance r to the boundary is held fixed and sets the scales. 
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FIG. 4. Illustration of the allowed values of u> (points) and the zeros of the argument of the delta- functions (line) in Eq. (pl[). 
The coordinates n and m are the summation indices. The line moves outward as u increases and a delta-function appears in 
the spectral density at values of u> where the line crosses a point with coefficients as given in Eq. (pl|). 
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FIG. 5. The predicted intensity J b s in arbitrary units as a function of to for boundary and bulk cases (i.e. for power-laws with 
abound = 1/2 and Qfbuik = 1/8, respectively). The corresponding three-dimensional case (a = 0) is also shown. Temperature 
(T = 50if) and finite resolution (A = 20meV) effects have been taken into account. 
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